Description of quantum dynamics of open systems based on collision-like models 
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Master equations in the Lindblad form describe evolution of open quantum systems that is com- 
pletely positive and simultaneously has a semigroup property. We analyze a possibility to derive 
this type of master equations from an intrinsically discrete dynamics that is modelled as a sequence 
of collisions between a given quantum system (a qubit) with particles that form the environment. 
In order to illustrate our approach we analyze in detail how a process of an exponential decay and a 
process of decoherence can be derived from a collision-like model in which particular collisions are 
described by SWAP and controUed-NOT interactions, respectively. 

I. INTRODUCTION 

The central issue in the study of evolution of open systems is the observed irreversibility of both classical as well 
as quantum dynamics. Open systems interacting with an environment gradually loose their information content and 
decohere, which means that after some time their states are to some extent determined by the initial parameters of 
the environment. Such behavior cannot be described by unitary transformations and has led to an introduction of a 
phenomenological dynamical postulate for open systems [1~3] - the semigroup property of the time evolution £t 

£t+s = £t£s ; Vi, s > , and hm £t = I , (1) 

where £t : 5(7i) S{T-C} is a transformation acting on the set of all possible states 5(7^) of a given quantum system 
corresponding to the Hilbert space 7Y. In this case it is possible that pure quantum states evolve into mixtures, and 
vice versa. 

The irreversibility is related to the non-existence of the inverse evolution £^^ . In particular, the inverse map can 
exist in the mathematical sense, but the resulting map does not describe a valid quantum evolution. The physical 
maps £t must satisfy several constraints [4]: They have to be linear, trace-preserving and they have to be completely 
positive for all t. 

The state of the system Qt is obtained by the application of the map on the initial state g, i.e. Qt = £t[Q\- Equations 
of motion that describe this type of dynamics has the so-called Lindblad form [1]: 

gt^G[Qt] = -t[H,Qt] + ^^c^i3{[Ac.,gt^i3] + [Kgt,^f3]) , (2) 

where Aq, = Aj^,, TrAo, = 0, TrA^A^ = daf3 for a,(3 = 1, . . . , — 1 and the coefficients Cap form a positive matrix. The 
Hamiltonian H is traceless, i.e. H = X^a^aAa. This differential equation is supplemented by the initial condition 

gt=o = g- 

Differentiating the time evolution Qt ~ £t [g] we derive a formal expression for the generator of the dynamics Q 

Qt = £t[g] = £t£t\gt] ^ Q = £t£t^ ■ (3) 

This expression in a strict mathematical sense has a meaning only when inverse transformations £^^ do exist for any 
time t. For general one-parametric set of completely positive maps £t the generator can be time-dependent, however, 
if these maps possess the semigroup property, then the generator is independent of time. 

II. DISCRETE DYNAMICAL SEMIGROUP 

The derivation of the particular master equation, which drives the evolution of open system, is usually based on 
the idea that the open system is a part of the larger closed system with the underlying unitary dynamics described by 
Schrodinger equation. However, after tracing out the environment, the induced set of maps £t essentially never fulfills 
the semigroup property. In order to obtain this feature of dynamics various approximations have to be applied [3,5], 
for instance the Born-Markov approximation, the weak-coupling limit, the mean-field limit, etc. Master equations 
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of the Lindblad form arc valid only under specific physical conditions, but still provide us with a very reasonable 
approximative picture of the exact dynamics of open quantum systems. 

In this paper we will present a novel method how to derive master equations. We will consider that the interaction 
between the system and its environment consists of bipartite collisions. For this we assume that the environment 
consists of N particles (quantum systems) of the same physical origin as the system under consideration. In our model 
we assume that the environment is initially "prepared" in a factorized state and the system interacts with each 
particle from the environment at most once. Each collision is described by a unitary transformation U, which induces 
a map £ (see Fig.l). As a result this collision-like model determines a discrete evolution described by powers of the 
map £, i.e. f " with n = 1, 2, . . . corresponding to the nth collision. The set of maps f„ := with So =1 obviously 
fulfills the semigroup property, i.e. 

SnSm = £n+m for all n, TO = 0, 1, 2, . . . (4) 

In other words the maps form a one-parametric discrete semigroup. 

The question is whether we may replace the discrete dynamics f „ with a continuous one £t such that £t^ = £n for 
times tn = TIT (r is some fixed time interval between subsequent collisions), and the maps £t satisfy the semigroup 
property. In other words the task is to interpolate a discrete set of points with a line in the abstract space of quantum 
maps. Obviously, on the time scale less than r the continuous evolution is different from the discrete one. However, 
our interest is to describe the overall evolution rather than a sequence of isolated collision. In this sense, the continuous 
dynamics is a good approximation of the discrete one. Our aim is to describe the intrinsically discrete collision model 
in the language of quantum dynamical semigroups, i.e. continuous master equations. 

We derive the master equation in the following way. Firstly, we will express the powers of the map f in a specific 
form. The discrete parameter n numbering the order of interaction of the system with environment particles will be 
replaced by a continuous parameter t = nr, i.e. n t/r. Correspondingly, we have a continuous set of maps £t. 
The question is, whether these maps are completely positive, and whether they form a semigroup. If, moreover, the 
inverse maps £f^ exist, then we can derive the master equation by using the expression for the generator in Eq. (3). 
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FIG. 1. A schematic visualization of a sequence of collisions between the system and particles from the environment. This 
interactions lead to a dynamics of the system that is described by a specific sequence of maps f„ that in the continuous limit 
can be described by a specific master equation. 




In what follows we will present two specific examples, for which the derived master equations correspond to an 
exponential decay, and to a decoherence process. It is important to understand which bipartite interactions under ly 
these type of processes. We will be interested only in two-dimensional systems, i.e. qubits. 



III. QUBIT FORMALISM 

In what follows we will use the so-called left-right formalism, or real matrix representation of quantum maps, and 
real vector representation of quantum states of a qubit. The state space of a qubit is a subset of four-dimensional 
linear space of hermitian operators. It follows that any density operator can be written as g = -\- f ■ a, where 
<^ — (fx, Cjf, Cz) are the Pauli operators. The positivity of g restricts the choices of vector r such that \ f\ < 1/2, i.e. 
the states of qubit form a three-dimensional sphere - Block sphere. 

The quantum evolution maps £ as well as the generators Q of quantum dynamics £ then correspond to 4x4 real 

matrices acting on the vectors (1/2, r). In particular, the evolution £ has the afHne form £ = ^ ^ ^ ^ > where T is 

some 3x3 matrix and t is the translation vector. Under the action of £ the Bloch sphere vectors f transform in the 
following way f^f'= Tf+t. The generator G of the dynamics has very similar form, only the first row of the matrix 
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vanishes, i.e. Q 







. The matrix elements are given by a relation [Q]jk = ^T^i'^jGWk]) for j, fc = 0, 1, 2, 3 and 



(To = /. 

For our purposes it will be useful to know how to rewrite this matrix form of the generator into the operator form 



0t 



1 ^ 

-i[H, gt] + - ^ Cjk {[(Jj,gt<Jk] + [(JjQt, tTfe]) , 



(5) 



j,k=i 



i.e. how to rewrite the coefficients [G]jk into the parameters hj and cjk- The hermitian matrix Cji~ can be rewritten as 

Cjk = djk — ie-jk, where djk = ■^{cjk+Ckj) is the real symmetric matrix and ejk = i^{cjk ^Ckj) is the real antisymmetric 
matrix. Using the operator expression of the generator (5), one can easily find the matrix 





Q = 



/ 

4623 - 2(i22 - 2^33 2di2 - 2h-i 

4631 2di2 + 2/l3 -2dii - 2d33 

V 4ei2 2^31 - 2/i2 2(^32 + 2hi 




(6) 



The inverse relations then read 
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and 



^ -fgl22-fgl33 + fg1ll 
. -[g]ll-[^]33 + [g]22 

"22 — 1 



d. 



-[g]ll-[g]22 + [g]33 
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di2 

d23 

di3 



[g]21-[g]l2 
[gl30 



ei2 



[g]l2 + [g]21 
[g]23 + [g]32 

4 

[g]l3 + [g]31 

4 
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It is important to note that these relations enable us to represent not only dynamical semigroups, but any time 
evolution of a qubit in the Lindbland-like form. In this general case the coefficients cjk will not be constant, but 
will explicitly depend on time. In this way we can derive more general master equations describing dynamics beyond 
Markovian approximation. 



IV. CASE STUDY I.: QUANTUM HOMOGENIZATION 



The quantum homogenization [6 8] is a process motivated by the thcrmodynamical process of thermalization. It 
describes a system-reservoir interaction in which the initial state of the system q is transformed into the state ^ 
determined by the state of the reservoir that is composed of N systems of the same physical origin as the system. 
The interaction between the system and the reservoir consists of bipartite collisions. Each collision is described by 
some unitary map U . In order to obtain a discrete semigroup describing the dynamics of the system, we assume that 
initially the reservoir is in a factorized state and that the system interacts with each system from the reservoir 
maximally once (see Fig.l). 

The homogenization approximate the evolution q®^^®'^ which is forbidden by the no-cloning theorem 

(see Ref. [9] and references therein). Let 6 be the parameter describing the quality of the homogenization in the 
following sense. After the homogenization process is complete, all systems are described by states, which belong to a 
^-vicinity of the state ^. Moreover, the homogenization requires the validity of the following relations corresponding 
to trivial homogenization Tri(C/^ ® ^J7^) = Ti2{U£,® £JJ'^) = ^. If we assume that the homogenization is independent 
on an initial state of the system qubit {q) and as well as on initial states of reservoir qubits (^) then for qubits U 
must posses the form of the partial swap operation (for more details see Ref. [6]) 

= cos rjl + i sin rjS (9) 

where S is the swap operation defined by the relation S\tp) (Si = ^ \tjj) for all \tp),\(j)). In what follows we will 
use the notation c = cos?7 and s = sinr]. 

Let us define the Bloch-sphere vectors for density operators, i.e. g ^ r and ^ <-> In terms of these vectors the 
partial-swap induces the map r^f'= c^r + s^t — 2cst x r, i.e. the superoperator £^ is represented by the matrix 
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Now we turn our attention on the derivation of the master equation that corresponds to a dynamics induced by 
a sequence of partial-swap interactions. By a suitable (unitary) substitution of the operator basis I , ax, cry, (Tz 
I, Sx, Sy, Sz, such that Sk — UaiJJ^ (UW — U^U — I), the matrix £^ can be rewritten into the form 




with A 



c 2sw 
-2sw c 



(11) 



where w is determined by the initial state of the reservoir, i.e. ^ ^ -^I + wSz- This change of the operator basis 
corresponds to a different choice of the x, y, z axes in the Bloch sphere. In what follows we will describe the process 
in this new basis, in which the powers of can be easily find 



cn 



1 

„n An ^ 

^ ^ 

2w{l - c^") c^" 



The powers of the matrix A can be found as follows: Using the identity cosarctanx = ^i-^-x'^ 
parameter co = arctan(2z/;s/c) we find 



(12) 

^ and defining the 



A = 



c 2sw 
—2sw c 



\/ + Aw^s'^ 



cos a; smw 
— sin uj cos uj 



and its powers are equal to 



(c2 + 4sW)"/2 



cos(ncj) sin(nw) 
— sin(7iti;) cos{nuj) 



(13) 
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The introduction of the continuous time is now straightforward. One has to replace n with t/r. The dynamics of 
open systems is usually characterized by two parameters: the decay rate Ti and the decoherence rate T2. In our case 
we can introduce these parameters as follows 



r.2t/T 



,-rit 



[c(c2+4sW)l/2]*/" 



ri = l/ri = -^lnc; 

= I/T2 = - i [In cVc2 +4s2u;2] , 



(15) 



and the continuous version of the homogenization process can be described by a one-parametric set of maps 



£t = 



^ 

e -^^cosilt e"^2*sinilt 
-e-^^^smQt e-^^*co8nt 



\ 2w{l 



-Fit 











-Fit 



(16) 



with frequency fl = uj/t describing the "rotating" part of the evolution (see Fig. 2). 
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FIG. 2. Visualization of the evolution of a system qubit due to its interaction with reservoir qubit in the process of quantum 
homogenization. The system is supposed to be initially prepared in a pure state (|0) + |l))/\/2. The qubits in reservoir are in 
prepared in the total mixture characterized by the values w = (left figure) and and in the ground state |0) characterized by 
w = 1/2 (right figure). 



The general rule is that before deriving of the master equation itself one has to check that the maps are valid 
quantum maps and that they fulfills the semigroup property. However, here we will directly derive the generator, 
and from its properties the character of transformation will become obvious. Using the method described in previous 
sections we find out that the derived generator is time independent and it takes the form 






-r2 -n 

-fl -Fa 

-Fi 



(17) 



Consequently, the master equation is of the Lindblad form with time-independent coefficients Cjk and hj. In particular, 
the non-vanishing parameters do read — Q/2, cn = C22 = ri/4, C33 = {2T2 — ri)/4, and C12 = —i2wri and the 
matrix of Cjk is positive whenever r2 > ^Fi, and \w\ < 1/2. The latter condition is trivially satisfied, because of 
the positivity of ^, i.e. \w\ < 1/2. The first condition is related to a more general result (for details see Ref. [10]), 
according to which the decay is maximally two times faster than the decoherence, i.e. Fi < 2F2, or equivalently 
T2 < 2Ti . Its validity in the case of homogenization can be checked by direct calculation. 

The positivity of the matrix cjk guarantees the complete positivity of the whole evolution. The master equation of 
the homogenization process reads 



Qt = -i^[Ss, Qt] + ^{SiqS^ + S2QS2 - 2g) + ^^{SsqSs - g) 
-iwTi{SigS2 - S2gSi + igSs + iSsg) , 



(18) 



where the operators S are given by Eq.(lO). Let us note that the parameter Q describes only the Hamiltonian part 
of evolution, which is in accordance with our expectations. The Hamiltonian part causes the rotation, whereas the 
other part of the generator induces a contraction into the fixed point ^. Only in the case when the homogenization is 
completed the map cannot be inverted, but strictly speaking this happens in the limit of time goes to infinity. 

Let us make the following choice of parameters: w = —1/2, Fi = 2F2 = 27. In this case the master equation takes 
a simple form 

Qt = i^[S3, g] + l[2S-gS+ - SS+g - gS-S+] , (19) 

where S± = ^{Si ± ^2). This is the well known equation describing the process of exponential decay with the 
Hamiltonian H = -^83 [4]. 



V. CASE STUDY II.: DECOHERENCE FROM COLLISIONS 



In this section we will analyze a collision- like dynamics which models a decoherence of a qubit. A more general 
treatment of decoherence will be presented elsewhere [11], while here we will concentrate to the derivation of the 
master equation. The task will be the same as before, except that instead of the partial-swap operation we will 
consider a partial CNOT (controUed-NOT operation), i.e. = cos 77/ + isinr^CNOT. The CNOT gate performs the ax 
rotation on the target qubit, when the control qubit is in the state |1). If the control is in the state |0), then the state 
of the target qubit is not changed. Unlike the swap operation, the controlled NOT is asymmetric under the exchange 
of qubits. Therefore, we have two different evolutions determined by the role of the system qubit: it can be either the 
target qubit, or the control qubit. 

Let us assume the system qubit acts as the target. In this case the interaction induces the map on the target qubit 
that reads 



^« = n n 1 -o„2.,, , (20) 
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where we used again the notation c — 00377,5 = sinry and stands for (1|^|1), i.e. it is the parameter that 
characterizes the initial state of reservoir qubits. For the calculation of the powers of this map we will use the same 
approach as before, because again we need to calculate the powers of 2x2 matrix A of the specific form 



A = 



-b a 



cos w sm uj 
— sin a; cosw 



= vV + 62X(a;) . 



(21) 



In the last equality we exploited the identity cosarctana; 



and we set w = arctan(6/a) to obtain the form 



suitable for evaluating the powers. Note that X(a;)" = X{niJ). In our particular case a = 1 — 25^^-^^ and 6 = 2cs^ii. 
That is, 



^5 



/ 1 

1 



V 



(l+4s2^n^oo)"/'^(nu;) 



(22) 



Following the same methods as in the previous section we can derive the generator of the process ^ in a very simple 
form 



/O 



-F 

VO -17 -F 



^ Qt = ~i-^[<^x, Qt\ + ■^{<yxQt<yx - Qt) ■ 



(23) 



As seen from its structure this generator describes a "pure" decoherence, i.e. a process of diagonalization of the state 
in the basis associated with the eigenvectors of the operator (j^- The parameters of the dynamics are defined in a 
similar way as before, i.e. = uj/t = arctan ^"^2^2^^^ and F = —7 In i^l — 'is^^^n^QQ. 

If we consider the system qubit to play the role of the control, then the dynamical map £^ is very similar to the 
previous case, i.e. 















1 

c2 + s2(^^,)^ C.S(1-K>«) 

-CS(1 - ((7^)e) C2 + S2(^^)^ 

1 



(24) 



where {<Jx)i = Tr^o-j,. Defining the parameters F 
write down the master equation 



- i In + s2(i7x){, Vt = ^ arctan 



cs(l— (o 



c2 + s2{o-a,> 



^ we can directly 



Qt = -'^-^[^z, Qt\ + -^{(^zQtf^z - Qt) ■ 



(25) 



Likewise in the previous case we obtain the dynamics describing the decoherence. Nevertheless, there is a difference, 
if the system qubit plays the role of the control, then the decoherence is observed in the basis associated with the 
eigenvectors of the cr^ operator, rather then in the basis of eigenvectors that was the case when the system qubit 
plays the role of the target. We see that irrespective on the fact whether the system qubit is a controlc or a target the 
sequence of partial CNOT collisions leads to decoherence processes. There are two difi'erences though. The first one is 
the basis in which the decoherence takes place and the second is the decoherence time. In particular, the decoherence 
rates depend on the initial state of the reservoir ^ in different ways. 



VI. CONCLUSION 



We have shown that using a simple collision-like model one can derive master equations for a qubit interacting 
with an environment. In our approach the discrete dynamics is described by a sequence of unitary transformations 
representing bi-partite interactions ("collisions"). As a result of the collision-like evolution the induced one-qubit 
dynamics is discrete, too. However, we have shown that in specific cases (the partial SWAP operation and the partial 
CNOT operation) this essentially discrete evolution can be substituted by a continuous dynamics. The resulting time 
evolution fulfills the semigroup property. Using the state-space description of the dynamics of a qubit we can say that 
the evolution of a qubit state that undergoes a sequence of collisions can be illustrated as a ordered sequence of points 
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in the Bloch sphere. We have shown how to connect these points with a smooth Une representing the continuous time 
evolution driven by a Lindblad master equation. 

Tlic presented approach can be generahzed to higher-dimensional quantum systems (qudits), as well as to collisions 
described by more general bi-partite unitary transformations U. The method essentially works for any interaction 
which induce an invertible map £ (i.e. det £ j^O). It turns out that semigroups £t = e~^* always contain only invertible 
maps, because for the (usually unphysical) map £-t the following identity holds £t£-t = e^^'e+^* = e" = X. 

Another (open) problem concerning the derivation of the master equation from a discrete collision-like dynamics 
is whether for all invertible mappings £, the semigroup generated by its powers can be always interpolated with a 
continuous semigroup of completely positive maps. This problem leads to the question about the structure of the 
set of quantum semigroups. It is worth to study which bi-partite interactions in the collision-like models could stand 
behind the known quantum processes described approximatively by quantum master equations. We believe that this 
approach of derivation of the master equations provides us with a new insight into the dynamics of open systems. 
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